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In this work we present a mathematical model for the mechanical response of the Brazilian Mario 
SCHENBERG gravitational wave (GW) detector to such waves. We found the physical parame- 
ters that are involved in this response assuming a linear elastic theory. Adopting this approach 
we determined the system's resonance frequencies for the case when six i-mode mechanical res- 
onators are coupled to the antenna surface according to the arrangement suggested by Johnson and 
Merkowitz: the truncated icosahedron configuration. This configuration presents special symme- 
tries that allow for the derivation of an analytic expression for the mode channels, which can be 
experimentally monitored and which are directly related to the tensorial components of the GW. 
Using this model we simulated how the system behaves under a gravitational sinewave quadrupolar 
force and found the relative amplitudes that result from this excitation. The mechanical resonators 
made the signal ~ 5340 times stronger. We found i -I- 1 degenerate triplets plus i non-degenerate 
system mode resonances within a band around 3.17 — 3.24kHz that are sensitive to signals higher 
than h ~ lO-^^Hz-^/^ when we considerate the effects of thermal noise only. 

Keywords: gravitational waves - instrumentation: detectors 



I. INTRODUCTION 

The existence of gravitational waves (GWs) was pre- 
dicted by Einstein applying his General Relativity The- 
ory to the vacuum in the weak field approximation 
0. GWs are local space-time curvature perturbations 
caused by accelerated masses. These perturbations travel 
through spacetime with the speed of light and can excite 
quadrupolar normal-modes of elastic bodies. 

Weber was the first one to propose feasible gravita- 
tional wave detectors, 43 years ago 0. However, so far 
nobody was able to convince the scientific community of 
any detection. The observable evidence of the existence 
of these waves comes from an indirect observation in the 
electromagnetic spectrum of a pulsar radio system 01 . 

In 1971 Forward suggested the use of a sphere as the 
antenna element of a resonant-mass detector Q . He ide- 
alized a sphere with nine sets of electromechanical strain 
transducers placed along great circle routes. The ten- 
sor gravitational radiation components would be deter- 
mined by five independent quantities from the nine trans- 
ducer outputs. Ashby and Dreitlein studied in detail the 
reception of GWs by an elastic self-gravitating spheri- 
cal antenna and Wagoner and Paik found the low- 
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est eigenvalues for the monopple and quadrupole modes 
of a uniform elastic sphere In the 1990's, Johnson 
and Merkowitz studied the antenna-transducer coupling 
problem and found an optimum configuration, which 
minimizes the number of transducers while keeping them 
in a symmetric distribution on the antenna: the trun- 
cated icosahedron (TI) configuration 7]. Magalhaes and 
collaborators showed that distributions with more res- 
onators also can have interesting symmetric properties 

In this work we extend Johnson and Merkowitz 's model 
to two-mode mechanical transducers. This kind of trans- 
ducer will be used on the Mario SCHENBERG spher- 
ical detector, which is being constructed in the Physics 
Institute of the Sao Paulo University with financial sup- 
port from the State of Sao Paulo Research Foundation 
(FAPESP). 

SCHENBERG has a 65cm-diameter and - 1150kg 
copper alloy [Cu(94%) Al(6%)] spherical antenna that will 
operate at temperatures below IK. These characteristics 
are shared with the detectors MiniGrail (The Nether- 
lands) and SFERA (Italy). Figure shows a schematic 
diagram of the Schenberg detector. 

In this work we consider an isolated sphere coupled 
only to mechanical resonators and predict the mechan- 
ical behavior of the system. We solve the displacement 
equations of the whole system and find an analytic ex- 
pression that relates the mode channels measurements to 
the GWs tensorial components. 
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FIG. 1: Schematic diagram of the Brazihan spherical GW 
detector Schenberg. 



II. 



THE SPHERE 



We will assume that the antenna has a perfect spher- 
ical symmetry and ignore the effects caused by the hole 
machined in the body to attach the suspension rod. We 
will show in subscction lV BI bv numerical simulations that 
this is a good approximation for the real case. 

Consider a spherical isotropic elastic body with radius 
R and density p in a field of forces. The strain suffered 
by a mass element dM is described by the displacement 
vector (where Xi represents the equilibrium 

position and x'^ denotes the position of dM after defor- 
mation). The strain tensor Uij is defined by 



{uij + Uj,i) , with i, j — 1, 2, 3. 



(1) 



There is a stress tensor associated to strain tensor by 
constitutive equations (generalizations of Hooke's law), 
given by 



(2) 



where A and /i are Lame coefficients, which describe elas- 
tic proprieties of a body (they are functions of the Pois- 
son's ratio and the Young modulus). By the application 
of the principle of conservation of linear momentum we 
find equilibrium equations given by 



where fi corresponds to the i-th component of an ex- 
ternal force density field and pui qq to an internal one. 
These equations are called Navier's equations and can be 
rewritten by application of equation |21 as 



(A + fi)uk,ki + P-Ui^jj + pUtflQ = /i, 



(4) 



The Navier's equations must be satisfied by a set of func- 
tions Ui = Ui{xi,X2TX^) which represent the displace- 
ments inside a region limited by a radius R. Therefore 
Equation 0] requires boundary conditions which will in- 
dicate that the surface is free to oscillate, defined as 



Ujaij =0, for 



(5) 



where rii = xi/r is the normal unit. 

We assume that it is possible to separate the spacial 
and temporal dependencies of the displacement vector at 
a position x in a time i, so we have 



Ui(x, t) 



,(i)*r„(x), 



(6) 



where Am{t) is the vibrational amplitude of the m-th 
normal mode and '4'm(x) is the eigenfunction of this 
mode. In fact, if the sphere is made of ordinary mat- 
ter {vsound ^ c) its radius is much smaller than the 
wavelength and its mechanical quality factor is high 



(Qr. 



i^niTm S> 1, where uJm and are the angular 



frequency and the energy time decay of m mode, respec- 
tively) such separation can be done. 

The eigenfunctions are normalized according to 



/ *„(x) • *,„(x)d3a; = iV„5„ 
Jv 



(7) 



where Nm is a arbitrary normalization factor for mode 
TO. If sphere is homogeneous (with constant p), we obtain 



in 



-i?^ V TO. 



(8) 



The equation of motion for the sphere as a forced har- 
monic oscillator can be obtained from Navier's equation 
and it is given by 



1 



pN„ 



[ *„,(x). Vf(x,t)d3 



(9) 



The eigenfunctions '4'm(x) lead information about the 
sphere's elastic proprieties (A and p), as well vibrational 
features (i.g. w^). 

Specific solutions for Equation|^can be obtained using 
spheroidal modes described as linear combinations of the 
quadrupolar spherical harmonics y„i(0, 0) to determine 
*m(x), which assumes the form lOj 



(3) 



*m(r,6i. 



[a[r)r + (3{r)RV]Y^ 



(10) 
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FIG. 2: Quadrupolar normal modes ol sphere. The darkest 
areas represent places of highest radial motion. 



where a(r) and are respectively the radial and tan- 
gential motion parameters. Figure |21 shows the distri- 
bution of the radial part of \J',„ on the sphere surface, 
the darkest regions representing the highest radial am- 
plitudes. Ashby and Dreitlein described air) and fi(r) 
function as Q 

3 R 
a{r) = piR-^j2iqr) +6p2—j2{kr) and (11) 
or r 



d 

P{r) = Pii2{qr) +Qp2-Q;^[rj2{kr)], 



{12) 



where j2 {x) is the spherical second order Bessel's function 
defined as |11|| 



3cosx 



sm X 



(13) 



and — pujq/{X + 2/i) and fc^ — pti^g/M respectively 
the squared transverse and longitudinal wave vectors. 
The value of the eigenfrequency ujq is determined from 
Equation which can be rewritten as 



Pi 



Pi 



dr 
12 





+ P2 
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2r2 



dr 



10- (fcr)2 
j2{qr) + 6p2 



1 d 
r dr 
d_ 

dr 



j2{kr)' 



0, 
= 0, 



where we consider only quadrupolar modes. The values 
of the normalization parameters pi and p2 are found from 
the normalization condition given by Equation 

Since we know the values of these parameters an ex- 
pression for effective gravitational force that acts on the 



TABLE I: Parameters for the SCHENBERG detector: intro- 
duced(l) and obtained from the present model(2). 



Description 


Symbol 


Value 


fl) 


sphere radius at 4K 


R 


0.3239m 


(1) 


sphere mass 


ms 


1147.85kg 


(1) 


sphere effective mass 




287.63kg 


(1) 


mass of resonator Ri 


ruR^ 


lO-^kg 


(1) 


mass of resonator R2 




0.054kg 


(1) 


density at 4K 


p 


8065.7kg/m^ 


(1) 


Young modulus at 4K 


E 


1.33 X 10"Pa 


(1) 


Poisson's ratio 


V 


0.364 


(1) 


linear contraction coefhcient^ 


^R/R 


334.52 X lO"'' 


(1) 


bulk Lame's coefficient at 4K 


fl 


4.8753 X lO^^Pa 


(1) 


shear Lame's coefficient at 4K 


X 


1.3049 X 10"Pa 


(1) 


sphere Q' 


Qs 


2.0 X lO'^ 


(1) 


first resonator Q 


Qri 


lO'' 


(1) 


second resonator Q 


Qr-2 


lO'' 


(2) 


degenerate mode frequencies 


fo 


3206.3Hz 


(2) 


normahzation parameter 1 


Pi 


-5.5654 


(2) 


normahzation parameter 2 


P2 


2.2758 


(2) 


radial parameter at r = _R 


aiR) 


2.8623 


(2) 


tangential parameter a.t r — R 


m) 


0.6598 


(2) 


Chi factor 


X 


0.6013 



"Tenderized mean (94%Cu 
''Mechanical quality factor 



sphere can be found as 



- 6%A1) from [T3. 



where 



X 



157r 



[pij2{qR) + 3p2j2{kR)] 



(14) 



(15) 



is a factor that determine the effective mass of the sphere, 
given by 



An 

12^ I Y 



R' 



4 ' 



(16) 



and hm = hm{h+{t),hy,(t),0^ip) denotes the spherical 
amplitude of gravitational wave that depends on the po- 
larization amplitudes and the orientation [0, relative to 
the lab coordinate system centered in the sphere's center 
of mass (SCM) 0. 

We solved the mathematical model for the sphere 
and found the specific physical parameters of the 
SCHENBERG detector, as is shown in Table U The 
physical dependencies between parameters can be seen 
in Figures |31 ^ and It is easy to notice that a{R) 
and x{R) have small dependence on the Poisson's ratio. 
Our tests shown that they are practically Young modulus 
independent. These features imply that little discrepan- 
cies on values used in our model will not affect hardly 
the results obtained in this work. 

Evidently, the gravitational force presented in Equa- 
tion^lis not the only force that acts on the sphere. The 
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FIG. 3: The normalization parameters dependencies on the 
Poisson's ratio i/. 
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FIG. 5: The radial a{R) and tangential P{R) displacement 
parameters and x(^) factor dependencies on v at the sphere 
surface {r — R). 
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FIG. 6: Unidimensional mechanical resonator with 2 modes 
attached to the sphere's surface. 



FIG. 4: The resonance frequency /o dependence on u and on 
the Young modulus E. 



III. THE SPHERE COUPLED TO 
RESONATORS 



resultant force has other components Fg that come from 
noise sources (e.g. Langevin's forces). These forces must 
be taken into consideration because they will contribute 
to the sum of the forces in Equation El Thus the sphere 
will be under the action of the forces F^m = + m- 
Many ways to minimize the noise forces contributions 
have been studied and a lot of what have been learned 
in the last decades with bar instruments has been reused 
and reorganized to be appHed to spherical detectors (e.g. 
vibration isolation systems, cryogenic devices, electro- 
magnetic isolation, etc). 

On the other hand, an isolated sphere is not a practical 
detector since the deformation it suffers due to the GW 
influence is too small to be detected. Something is needed 
to work as a impedance transformer. Thus we coupled 
secondary resonators to the antenna in order to amplify 
the signal. 



Suppose there are j resonators with 2 modes coupled 
to the sphere's surface, as it is shown in Figure |B| The 
resonator is assumed to displace only in the radial direc- 
tion. When the sphere's surface suffers a radial displace- 
ment described by the vector Zj under each resonator 
j, with masses mj^-^ and uir^ that will represent them 
henceforth, we have 

Zjit)^r,-J2Amm^{t), (17) 

m 

relative to the SCM, with fj being the unit radial vec- 
tor at the j-th resonator's position. Thus momentum 
is transferred to the j resonator, making mn-^ displace 
Iji + -^j from the SCM. Thus m/f^ moves a distance 
1j2 + Iji + from that point. 

The radial displacements of the sphere's surface un- 
der each resonator correspondent to each mode m can be 
collected in m vectors that describe the pattern of the 
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radial displacement for each mode. Those "pattern vec- 
tors" can be grouped in a "pattern matrix" B^^j , defined 

by la 



(18) 



where a = a{R) represents the radial displacement pa- 
rameter value at the sphere's surface. Using EauationllUI 
we can rewrite [T51 as 



Br, 



Y 

rr. 



(19) 



In other words, i?„i represents a spherical harmonic Ym 
behavior at Xj position Xj. 

The equations of motion of the system can be written 
in matrix form as 



msL g g 
mj^^aB^ mR-^I_ 

msL —'mR^aB_ 
"ij^,/ 





(20) 



where elements double underscored are matrices and 
single underscored ones are column vectors. We have 
B^ = [Bjm] as transpose of B_= [Bmj], L is the identity 
matrix and is the null matrix. The dimensions of these 
matrices depends only of the values of j and m. 

Equation 1201 was obtained assuming that the compo- 
nents of the system have high mechanical quality factor 
{Q > 10^). Consequently, the damping coefficients are 
negligible and the time decay of the oscillations tends to 
infinity. 



IV. THE RESONATORS' LOCATIONS 

Johnson and Merkowitz suggested that six resonator 
should be positioned in the center of the pentagonal faces 
of a truncated icosahedron (TI) concentric to the sphere 
0. The TI configuration has several special symmetries. 
The pattern matrix can be constructed in such a way 
that its pattern vectors will be orthogonal in pairs. It 
has the following properties (lOl 



BB ^ = 
Bl 



_3_ 
2^= 




(21) 
(22) 



and the following one, discovered during development of 
this work: 



= = 27r= 47r= 



(23) 



where 1 and \ are respectively a vector and a matrix in 
which all elements are the unity, and represents the null 
vector. 

These properties make the TI configuration a good op- 
tion to the resonator distribution on the sphere because it 
minimizes their effects on each other while using a small 
number of resonators, and make possible to obtain an 
analytic solution to the equations of motion, as shown in 
section Ivl 

Also, in TI configuration several resonators are affected 
by a specific sphere mode. So the sphere modes can be 
monitored by a linear combination of the six resonators 
outputs. Jonhson and Merkowitz have named these lin- 
ear combinations as mode channels, Qmit), and we have 
adapted their definition as 

9{t)^Ma2{t). (24) 

They are a direct amplified reading of sphere modes. 



V. A SOLUTION FOR THE EQUATION OF 
MOTION AND ITS ANALYTIC EXPRESSION 



Equation [5D| can be rewritten as 

2[r/w{t) + ujlxrfwit) = ZF(t). 



(25) 



where 



g 





w{t) = 




(26) 

and X, Y_ and Z_ are the original matrices from equa- 
tion 1201 following their original positions. In this way 
we solved the problem in a coordinate system based on 
center of mass, such 7 is the transformation matrix. 

We can normalize the equation 1251 bv multiplying it to 
(Xy)^^ through left and it became 



w{t) + LolMwit] ^ KF{t), 



(27) 



in which M = {^)^% and K = {Kir^Z. 

The next step is to obtain an orthogonal form of the 
equation of motion to separate the m + 2j harmonic oscil- 
lators from each other. That results in an equal number 
of linear equations and simplify the solution of the prob- 
lem. We assume that M is symmetric (and it really is) 
and therefore it can be diagonalized by using the relation 



UH'IU = D, 



(28) 



where f/ and D, represent respectively the matrix which 
columns are eigenvectors and the eigenvalues diagonal 
matrix of M. D_ must have positive and real values since 
frequencies must be real and positive. Thus we impose 
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that U_ is Hermitian {U}U_ = /) which assures M is pos- 
itively defined and D_ is positive and real. Therefore we 
can rewrite E7I as 

at)+ojimt)WKii(t), (29) 

where C — U}w. 

Now we can solve the equation of motion in the fre- 
quency domain where equation 1291 becomes 

{-10^1 + oJqR) C(^) = g^'KE.i'^)- (30) 

Since D. and / are diagonal matrices the matrix sum in 
the parentheses is invertible, implying that there is a ma- 
trix J (uj) that is equal to that sum and we have 

auj) = J(lo)U''K_F(u;). (31) 

To return to the original coordinates we reverse the trans- 
formations: 

qiicj) = jU^icj) = jUJ{u;)WKFiLj), (32) 
\ 92(w) / 

which represents the solution of the equation of motion 
in the frequency domain. 

A. Eigenvalues and eigenvectors 



B. The coupled mode frequencies 

Our tests for many i-mode identical resonators cou- 
pled to the sphere using TI configuration results 
that eigenfrequencies present i + 1 degenerated quin- 
tuplets and i non-degenerated modes. In Mario 
Schenberg's case we have three degenerate quintuplets 
(/- 3175.6, 3206.4, 3237.2Hz) and two isolated modes 
(/ - 3184.5, 3228.3Hz). Figure |VBl shows how resonance 
frequencies change as more transducers are coupled to 
the antenna. The number on the left corresponds to the 
number of resonators gradually added following the or- 
dering in Figure lVBI the number zero corresponds to the 
isolated sphere; the number one corresponds to resonator 
1 coupled to the sphere; the number two corresponds to 
resonators 1 and 2 coupled to the sphere and so forth. 
Dashed lines in the results with zero resonators indicate 
the frequencies measured of the real sphere at 4.2K. 

In Figure lVBI we present how the resonant frequencies 
should behave when we consider non-degenerate sphere 
modes, like those resultant from the measurements per- 
formed (the dashed lines we mentioned above). We con- 
cluded that there is a negligible effect from the suspen- 
sion hole on the symmetry proprieties of system so we 
decided to consider degenerate initial modes as a good 
approximation. 



Looking at eouation 1321 one can check that a good de- 
termination of eigenvalues and eigenvectors of M is es- 
sential to solve the problem. Moreover, the expressions 
must rigourously respect imposed conditions. From nu- 
merical results for Schenberg's case we noticed that U_ 
can be separated into two distinct groups: 




and 



gk^nk\ CkR^ I , for A: = 2,3,4, 
dkB^ 



(33) 
(34) 



where C/^ corresponds to three degenerate frequencies 
and U.i± to isolated frequencies, as is shown in subsection 
IV Bl The notation ± denotes eigenvectors which refer to 
up (-|-) and down (-) shifting of the frequency from the 
central degeneracy. The values of the constants ni±, di±, 
rife, Ck and dk are obtained from Eouation 1281 where 



D = 



Ai± 

XkL 



(35) 



and they are shown in appendix Thus we have 



C. An analytic expression to the solution of the 
equation of motion 

The special symmetries presented by the TI configu- 
ration enable us to find an analytic expression to the 
solution of the equation of motion |^. After a lengthy 
calculation we found an analytic expression for the mode 
channels, given by 



3(w) 



where 



1 4 

Erndk^kELsi'^) 




27r ^mR2 f- 



k=2 



3 a \ - , 

+ tT k , ly^Vkd^ 



b + a\ ~ N 
k\Ck ^ BFi 



k=2 
4 



277 hy/mn^ f-^^ \ h 
3 

27r h^/m^., 



-qkdk \ck + dkja + a ^)] RFj' , 

(38) 



9 = 



27r 



where 



jfj^ — Xi±ujq and ujI 



\ 2 



(36) 



(37) 



correspond to the resonant frequencies of the system. 



and 



b = 



ms 



7 



0) 
Q. 
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c 
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« 
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FIG. 7: Frequencies of tlie coupled modes. 




FIG. 8: Projected position of the resonators on the sphere's 
surface. 



The value of 77^ depends of Uk, Ck and dk values and it 
also is shown in appendix 1X1 



The above equation results from a combination of 
equations |^ and It represents the solution of the 
equation of motion and allows us to know the system's 
behavior under a gravitational wave excitation. 



VI. NUMERICAL SIMULATION OF THE 
SYSTEM 



In order to test the model we simulated a signal as a 
gravitational sinewave quadrupolar force with frequency 
equal to the degenerate normal mode frequencies of the 
sphere (~ 3206. 3IIz). The signal would be arriving at the 
local zenith so that its propagation direction coincided 
with the z axis of the reference system centered at SCM. 
Also, we assumed that the simulated gravitational wave 
had only "x" polarization ("+" polarization was null). 
Those aspects made only the second normal mode to be 
excited by the passing of the wave, as shown in Figure 
|51 We ignored any noise components in this simulation 
(Ff' = 0). 

Special symmetries presented by TI configuration al- 
low only some resonators to oscillate when only the sec- 
ond normal mode is excited, as shown in Figure ITUI The 
number on the right represent relative amplification to 
maximum amplitude of normal modes of the sphere. 

Richard showed (using the Principle of Energy Con- 
servation) that the ratio of displacements 0| 



(39) 



We apply this concept to Schenberg's case and found 

(40) 



\A\ 



5367, 
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0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 
t (msec) 

FIG. 9: Responses of the sphere's normal modes to the sim- 
ulated gravitational wave. Only the second normal mode was 
excited, as expected. 
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FIG. 11: Location of the transducers compared with second 
normal mode of the sphere. The darkest areas represent max- 
imum radial amplitude loci. 
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FIG. 10: Amplitudes of the second resonators {R2 
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FIG. 12: Mode channel outputs for the simulated signal. 



VII. AN ESTIMATE OF THE NOISE 
CONTRIBUTION ON THE MODE CHANNELS 



in resonant frequencies. This value is higher than that 
shown in Figure [TUI However as one can observe in Fig- 
urc llll no resonator is on a region of maximum amplitude 
which justifies the difference between the values. Fig- 
ure 111! presents a plane representation of second normal 
mode of sphere where darkest areas correspond to max- 
imum displacement loci. One can observe the resonator 
3 and 5 are placed on null radial movement areas. Res- 
onators 1 and 2 are placed close of maximum displace- 
ment points. An analogous analysis are done for all oth- 
ers normal modes and compatible results are obtained. 

From equation[21we obtain the mode channels outputs 
as showed in Figure [T^ 



We can rewrite Equation 1381 as 

gM = a^)F^ + ^l){^)Ffi + ^l)HFf2 (41) 

where ^(w) corresponds to the transfer function of the 

components on g^- Similarly, the Vl^^{uj) matrices 

represent the response functions of mode channel m to 

the noise source . 

... 

Assuming statistically independent noise sources and 
ignoring Tm on F^ , the mode channels spectral densities 
, S^{uj), can be written as 

Y^in^^icots^^'i, (42) 
j 
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FIG. 13: Spectral density of the contributions of the Langevin 
forces to the mode channels. 



FIG. 14: Sensitivity curve for the mode channels. 



with S^"' 



* and 5^^^ 



F^* FJI . As 

example one can calculate the noise contribution from 
Langevin forces, which have spectral densities well known 
and given by 



AkBTmiLj 



(43) 



where ks = 1.38 x 10~^'^J/K is Boltzman's constant, T 
is physical temperature of the system, and rrii and Qi 
are respectively the physical mass and the mechanical 
quality factor of the body i. 

From the calculation of the of the spectral density 
S!^ (lo), assuming no incident GW signal, we found that 
the contributions to each normal mode seem to have the 
same spectral densities, but with slightly different values, 
since each one of them receives different contributions 
from each of the resonators (See Figure IVIII) . These 
contributions are small compared to those originated 
from the sphere since the resonators' masses are much 
smaller than the sphere's mass {mg/Qs '^'^'ttiri/Qri ~ 
5.7x lQ^mR2lQR2, using the adopted values). The agree- 
ment among the values of 5*^^ {uj) is also a consequence 
of the fact that we have initially assumed that all the 
modes have the same transfer function, f {uj). The use 
of different transfer functions for each mode m may help 
the understanding on how asymmetries on the real sphere 
should affect the individual responses of the modes to 
GWs |1J]. 

The GW effective force (eq. I14f) that acts on the sphere 
(in the frequency domain) has spectral density 



The sensitivity curve for mode m is given through 



1 



(44) 



(45) 



Figure IVIII shows the sensitivity curves for the five 
quadrupolar normal modes of the sphere. We obtained 
this curve assuming that no signal was present and that 
only the usual noises in this kind of detector (Brownian 
noise, serial noise, back-action noise and electronic phase 
and amplitude noises) were disturbing the antenna. 

From reasonable estimates of the spectral densities of 
the noise sources it is possible to determine their con- 
tributions to the mode channels and, consequently, the 
values of hm (w), which allow for the location of the GW 
source in the sky, as well as the determination of the 
polarization amplitudes of the GW. 



VIII. CONCLUSION 

In this work we presented a model that describes the 
response to a gravitational wave of the spherical GW 
detector SCHENBERG coupled to six onc-dimensional 
two-mode transducers . 

The analytical expressions obtained as solutions of the 
equation of motion for the system are in good agreement 
with the numerical solutions. 

We found an expression fEauation l42|l which allows for 
estimates of the contributions of the noises to the mode 
channels. Such an expression is very much needed for a 
more realistic modelling of the SCHENBERG detector. 



APPENDIX A: EINGENVECTOR SOLUTIONS 

In order to determine the values of the constants 
Uk, Ck,dk,ni± and di± we imposed that C/ must be Her- 
mitian and 



MU = UD. 



(Al) 



10 



We separate the two groups of eigenvectors {U,ij. and f/^) 
and solve equation lAll using 

MRi± = Ui±Di±, (A2) 



12y3 + Ux^z - Sx^y^ + 81z^ - 54:xyz, (A13) 



(A3) 



where Di± = Ai± and D^. = are the eingenvalues 
of eigenvectors lLi± and respectively. So, by using 
the properties of the B matrix we obtain the following 
expressions for the constants: 



1 



di± = (a± + 4) , 



(A4) 



(A5) 



27r /a2 3 



27r 



^(.-|,-i 



' 9 5 ■ 



Those imply that 



Ai± = 1 — ad 



i± 



(A14) 



(A15) 



(A16) 



(A17) 



3 b f ^ ^ 2tt\ 



C2 



s§ + 4a;^ - 2xs3 - 2y 

6S3 



C3,4 



12S3 



+ 



S3 



where 



(A6) 



(A7) 



(A8) 



4(-l T iV3)x^ + (-1 ± iV^)4 - 4^X33 



(A9) 



(AlO) 



and 



Afe = 1 - —b Ck 

ZTT 



(A18) 



Although C3,4 {ck with fc = 3, 4) has imaginary terms 
its value is real, which satisfies Hermitian condition for 
eigenvectors. 

We also can define the quantities 



and 



rik 



(A19) 



(A20) 



These quantities help us to write an analytic expression 
for the mode channel outputs. 
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